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ABSTRACT 

We discuss our preliminary attempts to extend previous work on 2 x 2 Her- 
mitian octonionic matrices with non-real eigenvalues to the 3x3 case. 



1. INTRODUCTION 

In previous work [1,2; see also 3], we considered the real eigenvalue problem for 2x2 
and 3x3 Hermitian matrices over the octonions O. The 2x2 case corresponds closely to 
the standard, complex eigenvalue problem, since any 2x2 octonionic Hermitian matrix lies 
in a complex subalgebra C C O. The 3x3 case requires considerable care, resulting in some 
changes in the expected results. However, we also showed in [1] that there are octonionic 
Hermitian matrices which admit eigenvalues which are not real, and a complete treatment 
of the 2x2 case was given in [4]. Here, we discuss our preliminary results for the 3x3 case. 

Although we are able to obtain a 3rd-order characteristic equation for the (right) eigen- 
values in the 3x3 case, we have not been able to solve this equation, nor have we been able 
to extend our orthonormality results [1,2] from the real case. We therefore discuss several 
illustrative examples and make some conjectures regarding more general results. 

The 3x3 case is of particular interest mathematically because it corresponds to the 
exceptional Jordan algebra, also known as the Albert algebra. There have been numerous 
attempts to use this algebra to describe quantum physics, which was in fact Jordan's original 
motivation. More recently, Schray [5,6] has shown how to use the exceptional Jordan algebra 
to give an elegant description of the superparticle, which we have been attempting to extend 
to the superstring. Our dimensional reduction scheme extends naturally to this case [7] , and 
we believe it is the natural language to describe the fundamental particles of nature. 



1 Present address: Phoenix Corporate Research Laboratories, Motorola Inc., Tempe, AZ 85284, USA, 
r47569@email . sps . mot . com 
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Figure 1: The representation of the octonionic multiplication table using the 7-point projective 
plane, where we have used the conventional names {i,j,k,k£,j£,i£,£} for {e 2 ,...,e 8 }. 
Each of the 7 oriented lines gives a quaternionic triple. 

The paper is organized as follows. In Section 2 we briefly review the properties of 
octonions. We then consider 3x3 octonionic Hermitian matrices, deriving a characteristic 
equation for the eigenvalues in Section 3, and considering several examples in Section 4. 
Finally, in Section 5 we discuss our results. Some parts of this presentation have appeared 
in our previous work. 

2. OCTONIONS 

We summarize here only the essential properties of the octonions O. For a more detailed 
introduction, see [1] or [8,9]. 

The octonions O are the nonassociative, noncommutative, normed division algebra over 
the reals. In terms of a natural basis, an octonion a can be written 



The multiplication table is conveniently encoded in the 7-point projective plane, shown in 
Figure 1. The product of any two imaginary units is given by the third unit on the unique 
line connecting them, with the sign determined by the relative orientation. 

Octonionic conjugation is given by reversing the sign of the imaginary basis units 




(1) 



where the coefficients a q are real, and where the basis vectors satisfy e\ — 1 and 



e 2 q = -l (g = 2,...,8) 



(2) 



8 

a = a 1 ei — a q e, 

9=2 



(3) 
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Conjugation is an antiautomorphism, since it satisfies 

ab = ba 

The real and imaginary parts of an octonion a are given by 

Re(a) = i(a + a) Im(a) = ^(a - a) (4) 

The inner product on O is the one inherited from M 8 , namely 

a-b = J2 a q b q (5) 
i 



which can be rewritten as 



1 _ _ 1 _ 

a ■ b = -(ab + ba) = -(ba + ab) (6) 

2 2 



and which satisfies the identities 



a ■ (xb) = b ■ (xa) (7) 
(ax) ■ (bx) = \x\ 2 a-b (8) 



for any a,b,x G O. The norm of an octonion is just 



\a\ = v aa = -J a ■ a (9) 
which satisfies the defining property of a normed division algebra, namely 

\ab\ = \a\\b\ (10) 

The associator of three octonions is 

[a, b, c] = (ab)c — a(bc) (11) 

which is totally antisymmetric in its arguments, has no real part, and changes sign if any one 
of its arguments is replaced by its octonionic conjugate. Although the associator does not 
vanish in general, the octonions do satisfy a weak form of associativity known as alternativity, 
namely 

[b, a, a] = = [b, a, a] (12) 

The underlying reason for alternativity is Artin's Theorem [10,11], which states that that 
any two octonions lie in a quaternionic subalgebra of O, so that any product containing only 
two octonionic directions is associative. We will also have use for the associator identity 

[a, b, c]d + a[b, c, d] = [ab, c, d] — [a, be, d] + [a, b, cd] (13) 

for any a, b,c,d G O, which is proved by writing out all the terms. 

3. 3x3 OCTONIONIC HERMITIAN MATRICES 

In this section, we derive a characteristic equation for the (right) eigenvalues of A, which 
reduces to that of [1] for real eigenvalues. Unfortunately, we have been unable to solve this 
equation when the eigenvalues are not real, so that we have also been unable to investigate 
orthogonality and decomposition results analogous to those for real eigenvalues. We discuss 
this further in the next section, where we study several examples with intriguing properties. 
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a) Jordan matrices 

The 3x3 octonionic Hermitian matrices, henceforth referred to as Jordan matrices, form 
the exceptional Jordan algebra (also called the Albert algebra) under the Jordan product 2 

AoB:=\{AB + BA) (14) 

which is commutative, but not associative. A special case of this is 

A 2 = AoA (15) 

and we define 

A 3 :=A 2 oA = AoA 2 (16) 

Remarkably, with these definitions, Jordan matrices satisfy the usual characteristic equa- 
tion [13] 

A 3 - (trA)A 2 + a(A)A- {detA)I = (17) 

where cr{A) is defined by 

a(^):=i((tr^l) 2 -tr(^ 2 )) (18) 

and where the determinant of A is defined abstractly in terms of the Freudenthal product. 3 
Concretely, if 

(p a b\ 
a m c (19) 
b c n J 

with p, m, n G K and a, b, c G O then 

tr A = p + m + n 

o~(A) = pm + pn + mn — \a\ 2 — \b\ 2 — \c\ 2 (20) 



det A = pmn + b{ac) + b(ac) — n\a\ — m\b\ — p\c\ 

As shown originally by Ogievetsky [15], A, has 6, rather than 3, real eigenvalues, which 
furthermore fail to satisfy the characteristic equation (17). A complete, computer-assisted 
[2] treatment of this case was given in [1], and a somewhat more general analytic treat- 
ment was later given by Okubo [3]. As shown there, the eigenvalues naturally belong to 2 
distinct families, each containing 3 real eigenvalues. Furthermore, within each family, the 
corresponding eigenvectors lead to a decomposition of the form 

3 

A=Y,^(v a vl) (21) 
a=l 



2 The 2x2 octonionic Hermitian matrices form a special Jordan algebra since they are alternative [12]. 

3 The Freudenthal product of two Jordan matrices A and B is given by [14] 

A* B = Ao B - ^(Atr (B) + Btr (A)^j + ^ (tr (A) tr {B) - tr (A o B) 
The determinant can then be defined as 

det(A) = itr ((A*A)oA) 
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Furthermore, eigenvectors v a corresponding to different eigenvalues are automatically or- 
thogonal in the generalized sense 

(W) w = (22) 



b) Characteristic Equation 

Set 



Then the (right) eigenvalue problem 
becomes 



x 

y 1 (23) 



Av = vX (24) 



x (A — p) = ay + bz (25) 
y (A — m) = cz + ax (26) 
z (A - n) = bx + cy (27) 

As shown in [1], (24) admits solutions for which A is not real; further examples are given in 
Section 4. 4 Multiplying (25) on the right by (A — m) leads to 

(ay) (A - m) = x{\ - p) (A - m) - (bz) (A - m) (28) 

whereas multiplying (26) on the left by a 

a(y(X — m)) = a(cz) + \a\ 2 x (29) 

Subtracting these 2 equations immediately yields 

[a, y, A] = x ((A - p)(X - m) - |a| 2 ) - (bz)(X - m) - a{cz) (30) 

Similarly, multiplying (26) on the right by (A — p) and (25) on the left by a (or using 
symmetry) leads to 

[a, x, X] = y ((A — p)(X — m) — \a\ 2 ^j — (cz)(X — p) — a(bz) (31) 

We plan to multiply (30) by b on the left, (31) by c on the left, add, and then use (27). 
Before doing so, we first use (13) to write 

[b, x, (A — p)(X — m)] = [b, x(X — p), (A — m)] — [bx, (A — p), (A — m)} 

+ b [x, (A - p), (A - m)\ + [b, x, (A - p)\ (A - m) 
= [b,x(X-p),X] + [b,x,X](X-m) [ ' 

= [b, ay + bz, X] + [b, x, X] (A — m) 

as well as 

[b, bz, X] = [bb, z, X] + [b, b, zX] - b [l, z, X] - [b, b, z] X 
= -b[b,z,X] 



4 The quite different Jordan eigenvalue problem in the 3x3 case admits only real eigenvalues and was 
discussed in [7]. 
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(34) 



(35) 



Thus, returning to (30), we obtain 

fo [a, y, A] = (fox) ((A - p)(A - m) - \a\ 2 ) - \b\ 2 z(\ - m) - b(a(cz)) 

- [b, x, ( A — p) (A — m)] — b [fo, z, A] 
= (fox) ((A - p)(A - m) - |a| 2 ) - |5| 2 ^(A - m) - b(a{cz)) 

- [6, ay, A] - [6, x, A] (A - m) 

Similarly, (31) becomes 

c [a, x, A] = (cy) ((A — p)(X — m) — |a| 2 ) — \c\ 2 z(X — p) — c(a(bz)} 

- [c,ax,A] - [c, y, A] (A-p) 

Adding these last 2 equations, we obtain 

b [a, y, A] + c [a, x, A] = (fox + cy) (^(A — p){\ — m) — \a\ 2 ^j — \b\ 2 z(\ — m) — \c\ 2 z(\ — p) 

— b(a(cz)^J — c(li(bz)j 

- [c, ax, A] - [c, y, A] (A - p) - [fo, ay, A] - [fo, x, A] (A - m) 

(36) 

Finally, using (27), factoring out z, and rearranging terms leads to the generalized charac- 
teristic equation in the form 5 

z (A 3 - (tr^t) A 2 + a(A) A - detX) = b(a(cz)) + c(a(bz)) - (b(ac) + (ca)b)z 

+ fo [a, y, A] + [fo, ay, A] + [fo, x, A] (A - m) (37) 
+ c [a, x, A] + [c, ax, A] + [c, y, A] (A - p) 

If A is real, all the associators on the right- hand- side vanish, and we recover the gen- 
eralized characteristic equation given in [1]. The requirement in that case that the right- 
hand-side be a real multiple of z (since the left-hand-side is) then constrains z, resulting in 
precisely 2 values for that real multiple, and reducing (37) to 2 cubic equations, one for each 
family of real eigenvalues. 

While we find the form of (37) attractive, as there are no extraneous terms involving 
both z and A, we have so far been unable to further simplify (37) when A is not real. 

c) Alternate Approach 

We briefly describe another possible approach to finding the eigenvalues, which relies on 
the associator identity 

[f t ,f,A] := (v^v)X-v^vX) =0 (38) 

which follows for any octonionic vector v and A G O by alternativity, and which is further 
discussed in the Appendix of [4]. If v is a normalized right eigenvector of A with eigenvalue 
A, then 

v \Av) = v\v\) = (A)A = A (39) 



5 The extra terms multiplying z on the right-hand-side come from det A. 
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which yields an equation for A in terms of A and the components of v. A similar construction 
using the associator 

[v\ A, v] := {v ] A)v - v\Av) = (Av)K - v\Av) = (y\Av)f - v\Av) (40) 
leads for normalized eigenvectors to 

[v\ A, v] = (VM) f - v\v\) = ((v^v)X)f - (vK)\ = A - A = -2 Im(A) (41) 

Inserting (19) and (23) into (39) leads, after minor rearrangement using associators and 
(27), to 



\ x \ 2 + \y\ 2 ~ \ z \ 2 \ x \ 2 + \y\ 2 + \ z \ 2 ^ ^ 

which gives explicit expressions for the real and imaginary parts of A. The first term can be 
rewritten using cyclic permutations of {x, y, z) (and {a, c, b}), and the resulting expressions 
set equal to obtain 



x ■ (ay) + z ■ 


ibx) + p 


x 


2 




X 


2 



Re(A) = — \ X ( 43 ) 



and similar expressions obtained by cyclic permutation. 6 Finally, if v is normalized, the 
imaginary part of (42) reduces to 

Im(A) = [x, a, y] + [z, b, x] + [y, c, z] (44) 

We had hoped to use these various expressions to impose conditions on A which would in 
turn enable us to solve for A, but have not yet found a way to do so. 

4. EXAMPLES 

Without being able to solve (some version of) the characteristic equation in the 3x3 case, 
it is not possible in general to determine all the (non-real) eigenvalues of a given Hermitian 
octonionic matrix. It is therefore instructive to consider several explicit examples. 



a) Example 1 

Consider the matrix 



B 



( p iq kqs\ 

-iq p jq 
\-kqs -jq p J 



(45) 



where 



s = cos6> + k£sm9 (46) 

Note that B is quaternionic if 9 = 0. 

The real eigenvalues of B, and corresponding orthonormal bases of eigenvectors, were 
given in [2]. But B also admits eigenvectors with eigenvalues which are not real. For 



6 These expressions can also be obtained directly from (25)-(27) by taking the dot product with x, y, z, 
respectively, and using (7). 
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instance: 



A„ - p± qs ■ | ) S = 

J 

Xa=p±qs : v± = [ is | S± (47) 

A™ = V T : w± = \ —ks ] S± 
where 

St = { -f (48) 

These eigenvectors and eigenvalues reduce to the ones given in [2] when 6 — > 0. Somewhat 
surprisingly, these eigenvectors (when normalized) yield a decomposition of the form (21). 
Remarkably, they also yield a decomposition of the form 

3 

A= ]T KA Q ) ut (49) 

a=l 



We now describe some further properties of these eigenvectors. Each eigenspace with 
eigenvalues = p^f2qs is 1-dimensional (over R), so that the eigenvectors w± are essentially 
unique. By contrast, the eigenspaces with eigenvalues A^ = = p ± qs are 5-dimensional. 
Interestingly, though, for any given eigenvector such as u±, there is again an essentially 
unique eigenvector, in this case v±, which is orthogonal to it. Here "essentially unique" 
means unique up to a real multiplicative factor, and orthogonality can be defined either as 
v^w = or as (vv^)w = 0; these turn out to be equivalent in this case. 

There are also additional eigenvectors with eigenvalues of the form 

A = (p + p) - (3 U (50) 

where (3, p e M. These must satisfy 

32/3 2 = (^ V / 32p 2 - 7q 2 - 5q^j (llq - ^32p 2 - 7q 2 ^j (51) 
and therefore only exist provided that 

q 2 <p 2 < 4q 2 (52) 

In the case of equality, we recover the real eigenvalues given in [2]. For each admissible p 
and each of the two corresponding choices for f3, the eigenspace is 3-dimensional. We have 
not explored the properties of these eigenvectors in depth. 

All the eigenvalues discussed above have the form (50). While we suspect that there are 
no others, we have not been able to prove that this is the case. 
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b) Example 2 

A related example is given by the matrix 





( 


P 


qi 


§ks\ 


B = 




-qi 


V 






\ 




2J 


p J 



(53) 



with s again given by (46). We choose 9 such that 

3 3 



(54) 



B 



(55) 



resulting in 

V qi §(V5fe + 2Z) 

-qi p |j 

-§(V5* + 20 -|j p J 

The 2 families of real eigenvalues of B turn out to be {p± q,pT |(1 + V%),PT f (1 — ^r)}- 
Some eigenvectors for B corresponding to eigenvalues which are not real are 



V 



A 



"2 



V = (p - + y ^ 

A« 2 = (p - ^g) - 
A W1 = (p- ^g) - f Atf 

6 6 



A W2 = ip-^-q) + |h 



mi = 



■U2 



V\ = 



_ 3A; 

i + VEke 



V2 



W\ = 



W2 = 




(56) 



_ 3A; 
V5j -lit 

-7-VEke, 

( V5k + 2l \ 

3j 

\-3V5-3k£J 



However, we have been unable to find any decompositions of B involving these vectors. It 
is intriguing that, for instance, v\ is orthogonal to both u\ and w\ (in the sense of (22)), 
but that u\ and w\ are not orthogonal. In fact, we have shown using Mathematica that 
there is no eigenvector triple containing w\ which is orthogonal in the sense of (22). Unless 
w\ is special in some as yet to be determined sense, we are forced to conclude that neither 
(22) nor (21) are generally true for eigenvectors whose eigenvalues are not real. It is curious, 
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however, that the sum of the squares (outer products) of all six of these (normalized) vectors 
is indeed (twice) the identity! We will consider possible implications of this fact below. 

c) Example 3 

In all of the examples considered so far, the eigenvalues have been in the complex sub- 
algebra of O determined by the associator [a, b, c] (with a, b, c as in (19)). We now give an 
example for which this is not the case. 

Consider 

( V k -q(j-it-jt) 

C=\ -iq p q(l + k + l) | (57) 

\q(j -ii- j£) -g(l-fc-Z) p 
which admits an eigenvector 

' 3 







(58) 



with eigenvalue 

\ v =p + qlk (59) 

However, the associator takes the form 

[a, b, c 



= [i,(j-i£-j£),(l + k + l)]=2(l-k) (60) 



5. DISCUSSION 7 

As pointed out in [1], the orthonormality relation (22) is equivalent to assuming that 

vv^ + ... + ww^ = I (61) 

If we define a matrix U whose columns are just v, w, then this statement is equivalent to 

UU ] = I (62) 

Furthermore, the eigenvalue equation (24) can now be rewritten in the form 

AU = UD (63) 

where D is a diagonal matrix whose entries are the eigenvalues. Decompositions of the form 
(49) now take the form 

A = (UD)U ] (64) 
and multiplication of (63) on the right by W shows that 

(AU)U ] = (UD)U^ = A = A(Utf) (65) 

Thus, just as in [1], decompositions of the form (49) can be viewed as the assertion of 
associativity 

(AU)U ] = A(UU ] ) (66) 

We know of no way to express decompositions of the form (21) in similar language, which 
leads us to suspect that (49) is more fundamental. 



7 Much of this discussion also appeared in [4] . 
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This is somewhat less surprising when one realizes that 



((v\) u*) v = (vX) 



(67) 



for any (normalized) vector v and octonion A, due to an associator identity discussed in 
the Appendix of [4]. One can therefore construct matrices with arbitrary octonionic (right) 
eigenvalues, although such matrices will not be Hermitian. We conjecture, however, that the 
correct notion of orthogonality is 



not (22), to which it of course reduces if the eigenvalues are real, in which case it might be 
possible to decompose a matrix into non-Hermitian pieces of the form (67), and yet have the 
sum of these pieces be Hermitian. This is precisely what happens in the 2x2 case [4], as 
well as in Example 1. In any case, it is intriguing that this notion of orthogonality can be 
written as 



which explicitly involves A. 

Putting these ideas together, it would be natural to conjecture that all eigenvectors of 
a 3 x 3 octonionic Hermitian matrix come in families of 3, which form a decomposition in 
the sense that (66) is satisfied, and which are orthogonal in the sense of (68). However, 
this conjecture appears to fail for Example 2, as our claim that there is no orthonormal (in 
the sense (22)) eigenvector triple containing w\ rules out a decomposition (involving w\) in 
either the form (21) or the form (49). The reason for this is that the eigenvectors of B do 
not depend on the parameters p and q, which in turn means that the p and q parts of the 
decompositions can be treated separately. If a decomposition of either form were to exist, 
the terms involving p would then imply (62), contradicting our claim. (This is the reason 
that we did not reduce to the case p — 0, q — 1.) 

Nonetheless, for all known decompositions (62) does hold, and could be used to further 
simplify (66). One possible generalization would be for (62) to fail, but for a decomposition 
along the lines of (66) to hold. However, even (66) appears to fail for Example 2. 

There is, however, another intriguing possibility. Example 2 suggests that the eigenvec- 
tors of 3 x 3 octonionic Hermitian matrices may come in sets of 6 (or more), rather than 
in sets of 3. This would fit nicely with our recent result with Okubo [16] that, for real 
eigenvalues, it takes all 6 eigenvectors in order to decompose an arbitrary vector into a linear 
combination of eigenvectors, despite the fact that only 3 eigenvectors are needed to decom- 
pose the original matrix. Further evidence for this point of view is provided by the fact that 
the most general eigenvectors for the given eigenvalues of B have at most 2 free parameters, 
rather than the 4 degrees of freedom shown in [1] to exist for real eigenvalues, or the 8 for 
complex matrices. 

We therefore conjecture that, for any 3x3 octonionic Hermitian matrix, (66) should hold 
when suitably rewritten for a set consisting of n eigenvectors, where n presumably divides 
24, the number of (real) independent eigenvectors with real eigenvalues. However, we have so 
far been unable to write Example 2 in such a form. Whether or in what form orthogonality 
would hold in such a context is an interesting open question. 




(68) 




(69) 
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